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ABSTRACT 


This  report  clarifies  the  nature  of  the  dispersion  relation  for 
space-charge  waves  in  IMPATT  diodes.  It  is  demonstrated 
that,  in  the  usual  linear  approximation,  the  dispersion  rela¬ 
tion  is  always  cubic,  although  suitable  transformations  of  the 
basic  equation  appear  to  yield  a  quadratic.  The  implications 
of  this  point  are  discussed  in  regard  to  prior  results  and  for 
a  simple,  but  tractable,  generalization  of  these  results. 

In  addition,  possible  implications  for  the  TRAPATT-ARP 
controversy  concerning  the  explanation  of  anomalous  mode 
operation  of  avalanche  diodes  are  discussed. 


Accepted  for  the  Air  Force 

Franklin  C.  Hudson 

Chief,  Lincoln  I  aboratory  Office 
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DISPERSION  RELATIONS  FOR  IMP  ATT  DIODES 


This  report  is  concerned  with  a  clarification  of  the  basic  nature  of  IMPATT  wave  dispersion 
relations  which  apparently  are  not  widely  appreciated,  along  with  a  generalization  of  some  prior 
results.  In  the  first  seetion,  three  alternative  sets  of  basic  equations  will  be  derived  and 
diseussed. 

I.  BASIC  EQUATIONS 

One  basie  set  of  equations  for  one-dimensional  interactions  is  the  continuity  equations  for 

2 

electron  and  hole  currents  along  with  the  Poisson  equation.  These  are 
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where  J  =  —  qn  V  and  J 
n  ^  n 1  p 


the  electron  and  hole  densities,  q  is  the  electronic  charge, 


qp|Yp|  are  the  electron  and  hole  current  densities,  n  and  p  a  re 

Y  and  \  are  the  magnitude's  of 
n  p 

the  electron  and  hole  velocities,  rv  and  (3  are  the  electron  and  hole  ionization  rates,  N  and  \ 

o  a 

are  the  time-invariant  donor  and  acceptor  densities  as  determined  by  the  material  doping,  <  is 
the  permittivity,  and  E  is  the  electric  field.  This  set  involves  three  equations  and  three  un¬ 
knowns  (n.  p,  K  ),  and  hence  the  dispersion  relation  can  be  expeeted  to  be  cubic  when  the  equa¬ 
tions  arc  linearized. 

A  second  set  ean  be  derived  in  at  least  two  ways.  In  the  first  way,  Eq.  (2)  is  subtracted 
from  Eq.  (1)  to  yield 
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If  Eq.  (3)  is  differentiated  partially  with  respect  to  time,  the  result  is 
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Substituting  Eq.(5)  into  Eq.  (4),  we  obtain 
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which  may  be  partiallv  integrated  with  respect  to  x  to  yield 

BE. 
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where  J^,  called  the  total  eurrent  density,  represents  the  sum  of  the  particle  and  displacement 
current  densities  and  is  solely  a  function  of  time.  Equations  (1),  (2),  and  (7)  constitute  the  al¬ 
ternative  set  of  three  equations  in  three  unknowns  (n,  p,  E^)  for  which  a  cubic  dispersion  relation 
is  expeeted  when  the  equations  are  linearized.  The  second  method  of  derivation  recalls  that  t  he  * 
Maxwell  e qua t ions, 
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must  be  applicable  to  E  in  Eq.  (3).  Taking  the  divergence  of  Eq.  (8),  we  obtain 
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because  V  •  V  X  H  -  0  is  a  vector  identity.  For  one- dimensional  interactions  (V  x  E 

a  -  a  -  a  ,  ° 

E  x  ,  J  =  J  x  ,  J  =  J  x  ),  Eq.  (10)  becomes  Eq.  (6)  and  hence  Eq.  (7). 
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In  the  small- signal  time-harmonic  case  where 
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the  equations  for  the  time-harmonic  components  are 
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(where  g  =  cv  n  |  V  +  £  p  |  V  1  f  a'E.N  |  V  |  +  p'E.P  |V  | }  or  Eqs.  (15),  (16)  and 
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(where  OJ^  /fix  =  0).  It  is  important  to  note  that  Eq.  (17)  can  be  derived  from  Eqs.  (15),  (16), 
and  (18).  This  could  not  be  done  in  the  large-signal  case  because  Eqs.  (1),  (2),  and  (6)  combine 
to  give 
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where  F(x)  is  an  unknown  function  [unknown  because  information  not  contained  in  Eqs.  (1),  (2), 
and  (6)  must  be  invoked  to  prove  that  F(x)  =  (q/e)  (Nq  —  N^)]. 


2 


Some  investigators  regard  the  derivation  of  Eq.  (17)  from  Eqs.  (15),  (16),  and  (18)  as  indica 
tion  that  there  are  only  two  independent  equations  [Eqs.  (15)  and  ( 1  6) j  and  hence  that  the  disper- 
sion  relation  must  be  quadratic.  They  apparently  neglect  to  include  Eq.  (18)  in  their  equation 
count.  It  is  the  form  of  the  third  alternative  set  of  equations,  which  we  now  derive,  which  ap¬ 
parently  encourages  the  preceding  mistaken  belief. 

If  Eq.  (17)  is  differentiated  with  respect  to  x  to  obtain 
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and  Eqs.  (15)  and  (16)  are  substituted  in  Eq.  (19),  the  result  is 
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In  the  very  special  case  where 
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are  constants  if  E  and  hence  o  ,  3  ,  o’  and  3'  are  independent  of  x  as  in  the  case  of  a  PIN 
configuration.  Because  Kq.(21)  is  a  second-order,  inhomogeneous  differential  equation  for  Et 
some  investigators^  take  this  as  further  evidence  that  the  dispersion  relation  must  be  quadratic. 
This  view  overlooks  the  point  that  the  solution  to  Eq.  (21)  will  include  a  constant  term  (i.e.,  in¬ 
dependent  of  x)  due  to  the  presence  of  the  constant  "forcing  function,"  B.l  ^ ,  in  the  differential 
equation.  This  constant  term  in  the  solution  corresponds  to  the  K  0  root  (when  substituted  into 

-  j 

e  *  )  of  the  cubic  dispersion  relation.  If  the  cubic  dispersion  relation  did  not  possess  a  K  0 

root,  the  reduction  of  the  system  of  three  equations  [Eqs.  (15),  (16),  and  (17)]  to  the  single  second- 
order  differential  equations,  Eq.  (21),  would  not  be  possible.  Thus,  the  solutions  to  the  three 
coupled  first-order  equations  have  the  form 
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which,  if  one  uses  the  information  that  K,  -  0,  becomes 

-jK-x 

{n1,p1,E1}  =  Ai  +  A2e  +  A3e  ,  (23) 

which  is  the  form  of  the  solution  for  Eq.  (21).  The  basic  third- order,  homogeneous  differential 
equation  is  derived  in  the  Appendix. 


II.  DISPERSION  RELATION 

Manas se  and  Shapiro^  have  explored,  at  length,  a  generalization  of  Misawa's  dispersion 
2 

relation  which  is  more  physically  realistic,  while  still  remaining  tractable,  because  it  takes 
into  account  the  difference  in  electron  and  hole  parameters.  This  dispersion  relation  is^ 
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However,  Manasse  and  Shapiro  ignore  the  double  root  at  K  0,  of  Eq.(24),  and  do  not  regard  it 

as  part  of  the  true  solution.  Rather  they,  as  well  as  others,  considered  it  to  be  quadratic,  ap- 

1 

parentlv  for  some  of  the  reasons  discussed  in  the  first  section. 

We  consider  now,  briefly,  a  simple  and  tractable  generalization  of  the  preceding  results. 

We  note,  from  Sze  and  Cibbons,^  that  n  *  /a  need  not  equal  il'  /li  ,  and  in  fac  t  differ  bv  a  fac  tor 

o  o  o  o 

of  approximately  two  for  silicon.  Thus,  from  Eqs.(15),  (46),  and  (47),  we  find  for  traveling- 
wave  solutions 
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where  V  =  |V  |\  j/d.  When  a  /(l  -  c v  •  //>*  ,  Eq.  (26)  reduces  to  Eq. (24)  except  for  an  extra 
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factor  oi  d  in  the  third  term,  .  .  .  r  jd(rv  +  G  d )  1 ,  of  the  left-hand  side  of  Hq.  (26)  whose  absence 

J  o  o  ^ 

in  Eq.  (24)  is  presumably  due  to  a  typographical  error.  The  solutions  to  Eq.  (26)  are 
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The  inverse  dispersion  relation,  cc  =  g(K)  is  obtained  from 
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A  deseription  of  o  '  .1  +  /i’J  ,  given  ,1  ,  ,  is  required  to  make  use  of  the  above.  This  require 

r  o  no  o  po  dc.  n  n 

a  straightforward  generalization  of  the  argument  given  in  Appendix  II  of  Ref.  3.  The  result  is 

from  an  averaging  procedure  for  the  de  results  and  yields  for  o'.)  +  B' .]  =  G, 
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where  e  =  pQ/cy  ,  e  =  fi'0/a  '0>  and  L  is  the  length  of  the  avalanche  region.  When  e  -  c,  the  result 

reduces  to  Manas se  and  Shapiro’s  .1  =  J  ,  /a  L. 
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The  preceding  result  was  obtained  by  using  o'.T  +  11  =  rv 1  [ «T  .  +  (e  —  1)  J  1  (recall  that 
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III.  CONCLUSIONS 

We  c  onsider  now  the  application  of  some  of  the  preceding  results  to  a  controversial  topic. 

It  has  been  demonstrated  that  IMPATT  diode  operation  is  described  by  a  third-orcler  homoge¬ 
neous,  differential  equation  which  gives  rise  to  a  cubic  dispersion  relation.  This  conclusion  is 
contrary  to  the  widely  held  belief,  by  workers  in  the  field,  that  the  basic  equation  and  resulting 
dispersion  relation  are  of  second  order.  This  distinction  does  not  appear  to  he  of  practical  im¬ 
portance  for  small- signal  theory  of  IMPATT  diodes,  since  all  three  roots  of  the  dispersion  re¬ 
lation  (0,  +k  ,  —  k  )  are  used  in  practice  bv  everyone  to  calculate  diode  properties, 
m  m  - — - 1 - - - 

There  may  be  important  consequences  in  large- signal  ayalanche  diode  theory.  At  the  pres¬ 
ent  time,  proponents  of  the  TRAP  ATT  mode  theory^  (as  an  explanation  of  anomalous  mode  ava¬ 
lanche'  diode  behavior)  argue  that  they  have  solved  the  general  system  of  equations  by  computer 
simulation  and  see  only  TRAPATT  mode  operation.  Other  workers^  claim  that  there  exists  an 
additional  high-efficiency  mode,  avalanche- resonance  pumped  (ARP),  which  they  obtain  ex¬ 
perimentally  and  which  docs  not  hat  e  the  waveforms  (of  J  and  voltage  vs  time)  predicted  by 
TRAPATT  theory. 

The  solution  to  the  controversy  may  lie  in  the  observation  that  the  computer  simulations  of 
TRAPATT  mode  operation  assume  given  waveforms  for  .1^  (typically  a  step- modulated  sine  wave) 
which  is  equivalent  to  solving  the  system  as  though  it  were  described  completely  by  a  second- 
order  partial-differential  equation.  It  seems  plausible  to  this  author  that  if  were  not  specified, 
but  treated  as  the  unknown  it  truly  is,  a  more  general  solution  might  be  obtained  which  would 
contain  ARP  and  TRAPATT  modes  as  separate  possible  solutions.  It  maybe  argued  in  opposition 


5 


to  this  view  that  it  is  known  that  .1^,  in  the  one-dimensional  case,  is  solely  a  function  of  time 
and  that  the  assumed  time  dependence  for  can  be  arrived  at  by  reasonable  physical  arguments. 
This  is  true  but  not  conclusive,  since  it  is  pointed  out  on  page  2  that  in  the  large-signal  case 
Eqs.  (1),  (2),  and  (6)  are  not  equivalent  to  (in  the  sense  that  they  cannot  be  used  to  derive)  Eqs.  (1), 
(2),  and  (3).  Thus,  assumptions  on  the  form  of  d  (t)  may  omit  other  physically  possible  but  not 
obvious  solutions. 
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APPENDIX 

DERIVATION  OF  BASIC  EQUATION 


The  basie  set  of  three  linearized,  first-order,  eoupled  differential  equations,  Eqs,  (15), 
(16),  and  (17),  combine  to  yield  a  third-order,  homogeneous  differential  equation  which  we  dis¬ 
play  in  this  appendix. 

Eqs.  (15)  and  (16)  may  be  written  in  the  form 
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I  rom  Eqs.  (34)  and  (3  5)  it  can  be  shown  that 
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Thus,  if  Eq.  (17),  which  we  repeat  here. 
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